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We derive the mass-radius relation of relativistic white dwarf stars (modeled as a self-gravitating 
degenerate Fermi gas at T = 0) in a D-dimensional universe and study the influence of the dimension 
of space on the laws of physics when we combine quantum mechanics, special relativity and gravity. 
We exhibit characteristic dimensions D = 1, D = 2, D = 3, D = (3+VT7)/2, D - 4, D = 2(l + v / 2) 
and show that quantum mechanics cannot balance gravitational collapse for D > 4. This is similar 
to a result found by Ehrenfest (1917) at the atomic level for Coulomb forces (in Bohr's model) 
and for the Kepler problem. This makes the dimension of our universe D = 3 very particular with 
possible implications regarding the anthropic principle. We discuss some historic aspects concerning 
the discovery of the Chandrasekhar (1931) limiting mass in relation to previous investigations by 
Anderson (1929) and Stoner (1930). We also propose different derivations of the stability limits of 
polytropic distributions and consider their application to classical and relativistic white dwarf stars. 

PACS numbers: 05.90.+m; 05.70.-a; 95.30.-k; 95.30.Sf 



INTRODUCTION: SOME HISTORIC 
ELEMENTS 



The study of stars is one of the most fascinating topics 
in astrophysics because it involves many areas of physics: 
gravitation, thermodynamics, hydrodynamics, statistical 
mechanics, quantum mechanics, relativity... and it has 
furthermore a very interesting history. In his classical 
monograph The Internal Constitution of the Stars, Ed- 
dington (1926) [l[ lays down the foundations of the sub- 
ject and describes in detail the basic processes that gov- 
ern the structure of ordinary stars. At that stage, only 
elements of "classical" physics are used and difficulties 
with such theories to account for the structure of high 
density stars such as the companion of Sirius are pointed 
out. Soon after the discovery of the quantum statistics by 
Fermi (1926) @| and Dirac (1926) 0, Fowler (1926) [j 
uses this "new thermodynamics" to explain the puzzling 
nature of white dwarf stars. He understands that low 
mass white dwarf stars owe their stability to the quantum 
pressure of the degenerate electron gas 52] . The resulting 
structure is equivalent to a polytrope of index n = 3/2 
so that the mass-radius relation of classical white dwarf 
stars behaves like MR 3 ~ 1 (Chandrasekhar 1931a @]). 
The next step was made by Chandrasekhar, aged only 
nineteen, who was accepted by the University of Cam- 
bridge to work with Fowler. In the boat that took him 
from Madras to Southampton [9J] , Chandrasekhar under- 
stands that relativistic effects are important in massive 
white dwarf stars and that Einstein kinematic must be 
introduced in the problem. In his first treatment (Chan- 
drasekhar 1931b [l(|), he considers the ultra-relativistic 
limit and shows that the resulting structure is equivalent 
to a polytrope of index n = 3. Applying the theory of 
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polytropic gas spheres (Emden 1907 [ll|), this leads to 
a unique value of the mass M c that he interprets as a 
limiting mass, nowdays called the Chandrasekhar limit 
[53j. The complete mass-radius relation of relativistic 
white dwarf stars was given later (Chandrasekhar 1935 
|13l |) and departs from Fowler's sequence as we approach 
the limiting mass. At this critical mass, the radius of 
the configuration vanishes. Above the critical mass, the 
equation of state of the relativistic degenerate Fermi gas 
of electrons is not able to balance gravitational forces 
and, when considering the final evolution of such a star, 
Chandrasekhar (1934) [1J| "left speculating on other pos- 
sibility" . Chandrasekhar's result was severly criticized by 
Eddington (1935) [15j who viewed this result as a reduc- 
tio ad absurdum of the relativistic formula and considered 
the combination of special relativity and non-relativistic 
quantum theory as an "unholly alliance" . Because of 
the querelle with Eddington, it took some time to real- 
ize the physical implication of Chandrasekhar's results. 
However, progressively, his early investigations on white 
dwarf stars were extended in general relativity to the case 
of neutron stars (Oppenheimer & Volkoff 1939 [16|]) and 
finally led to the concept of "black holes" , a term coined 
by Wheeler in 1967. Without the dispute with Edding- 
ton, this ultimate stage of matter resulting from gravi- 
tational collapse could have been predicted much earlier 
from the discovery of Chandrasekhar [17| ■ 

Although these historic elements are well-known, it is 
less well-known that the concept of a maximum mass 
for relativistic white dwarf stars had been introduced 
earlier by Anderson (1929) and Stoner (1930) 
[54J . These studies are mentioned in the early works of 
Chandrasekhar but they have been progressively forgot- 
ten and are rarely quoted in classical textbooks of as- 
trophysics. These authors investigated the equation of 
state of a relativistic degenerate Fermi gas and predicted 
an upper limit for the mass of white dwarf stars. Stoner 
(1930) [19( uses a uniform mass density to model the 



star while Chandrasekhar (1931) [KJ considers a more 
realistic n = 3 polytropc. However, as noted in Chan- 
drasekhar (1931) [l0(, the value of the limiting mass 
found by Stoner with his simplified model is relatively 
close to that obtained with the improved treatment. In- 
terestingly, Nauenberg (1972) [21j introduced long after 
a simplified treatment of relativistic white dwarf stars in 
order to obtain an analytical approximation of the mass- 
radius relation. It turns out that this model, which gives 
a very good agreement with Chandrasekhar's numerical 
results (up to some normalization factors), is similar to 
that introduced by Stoner. 

Leaving aside these interesting historical remarks, the 
result of Chandrasekhar [10| concerning the existence of a 
limiting mass is very profound because this mass can be 
expressed in terms of fundamental constants, similarly 
to the Bohr radius of the hydrogen atom [22|]. Hence, 
the mass of stars is determined typically by the following 
combination 
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where G is the constant of gravity, h the Planck constant, 
c the velocity of light and H the mass of the hydrogen 
atom (Mq is the solar mass). This formula results from 
the combination of quantum mechanics (h), special rel- 
ativity (c) and gravity (G). Since dimensional analysis 
plays a fundamental role in physics, it is of interest to 
investigate how the preceding results depend on the di- 
mension of space D of the universe. In our previous inves- 
tigations 23, 24, 25], we considered the case of classical 
white dwarf stars in D dimensions and found that they 
become unstable in a space of dimension D > 4. In that 
case, the star either evaporates or collapses. Therefore, 
quantum mechanics cannot stabilize matter against grav- 
itational collapse in D > 4 contrary to what happens in 
D = 3 @, 0, Q ■ Interestingly, this is similar to a result 
found by Ehrenfest [26] at the atomic level for Coulomb 
forces in Bohr's model and for the planetary motion (Ke- 
pler problem). The object of this paper is to extend these 
results to the case of relativistic white dwarf stars and ex- 
hibit particular dimensions of space which play a special 
role in the problem when we combine Newtonian gravity, 
quantum mechanics and special relativity. We shall see 
that the problem is very rich and interesting in its own 
right. It shows to which extent the dimension D = 3 
of our universe is particular, with possible implications 
regarding the anthropic principle [23]. We note that a 
similar problem has been considered in [28( on the basis 
of dimensional analysis. Our approach is more precise 
since we generalize the exact mathematical treatment of 
Chandrasekhar 13] to a space of dimension D. A con- 
nection with other works investigating the role played by 
the dimension of space on the laws of physics is made 
in the Conclusion. In Appendix IBl we propose different 
derivations of the stability limits of polytropic spheres 
and consider applications of these results to classical and 
relativistic white dwarf stars. 



II. THE EQUATION OF STATE 

Following Chandrasekhar [13j , we model a white dwarf 
star as a degenerate gas sphere in hydrostatic equilib- 
rium. The pressure is due to the quantum properties of 
the electrons and the density to the protons. In the com- 
pletely degenerate limit, the electrons have momenta less 
than a threshold value po (Fermi momentum) and their 
distribution function is / = 2jhP where h is the Planck 
constant. There can only be two electrons in a phase 
space element of size hP on account of the Pauli exclu- 
sion principle. Therefore, the number of electrons per 
unit volume is 
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where Sd = 2n D / 2 /T(D/2) is the surface of a unit sphere 
in _D-dimensions. The mean kinetic energy per electron 
is given by 



* = \J fe(p)dp = jg§ J e(p)p"-'dp, (3) 



where e(p) is the energy of an electron with impulse p. 
In relativistic mechanics, 
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The pressure of the electrons is 

lf.de 2S D f° D de 

p= dJ fp d-p dp = Dh^J Q p T P dp - (5) 

Using Eq. ((4|), the pressure can be rewritten 
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Finally, the mass density of the star is 

P = nfiH, (7) 

where H is the mass of the proton and /i the molecular 
weight. If we consider a pure gas of fermions (like, e.g., 
massive neutrinos in dark matter models), we just have 
to replace fj,H by their mass m. 

Introducing the notation x — py/mc, we can write the 
density of state parametrically as follows 
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The function f(x) has the asymptotic behaviors 



The density can be expressed as 
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The classical limit corresponds to i < 1 and the ultra- 
relativistic limit to a; > 1. Explicit expressions of the 
function f(x) are given in Appendix [X] for different di- 
mensions of space. 



III. THE CHANDRASEKHAR EQUATION 

For a spherically symmetric distribution of matter, the 
equations of hydrostatic equilibrum are 
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They can be combined to give 
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Expressing p and P in terms of x and setting y 2 = 1 + x 2 , 
we obtain 
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We denote by Xq and j/o the values of x and y at the 
center. Furthermore, we define 



r = a?], y = y <f>, (17) 

yl = l + xl (18) 



SA 2 \ 1/2 I 



SnG 



(D-l)/2> 



D^J By', 

Note that the scale of length a is independent on yo for 
D = 1. Substituting these transformations in Eq. ((THJ), 
we obtain the Z?-dimensional generalization of Chan- 
drasekhar's differential equation 
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with the boundary conditions 

0(0) = 1, <j>'(0) = 
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The radius R of the star is such that p(R) = 0. This 
yields 
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where the central density 

p = Bx° = B(y 2 \) D ' 2 . (23) 

Finally, we find that the mass is related to yo by 
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Note that yo does not explicitly enter in this expression 
for D = 3 but it is of course present implicitly. 



IV. THE CLASSICAL LIMIT 

In the classical case i < 1, we find that the equation 
of state takes the form 
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Therefore a classical white dwarf star is equivalent to a 
polytrope of index [23| [55| : 
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Polytropic stars are described by the Lane-Emden equa- 
tion [11J| . This can be recovered as a limit of the Chan- 
drasekhar equation. Forx <C 1, we have yo — 1+^Xq. We 
define 6 = 2 -l/y§. To leading order, <\> = l-(x|-0)/2. 
Setting ^ = i/2r] and combining the foregoing results, we 
find that Eq. (|19|) reduces to the Lane-Emden equation 
with index D/2: 
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Note that the condition at the origin is 9(0) = Xq instead 
of 9(0) = 1 as in the ordinary Lane-Emden equation. 
However, using the homology theorem for polytropic 
spheres [29]], we can easily relate 9 to 9 D / 2 i the solution 
of the Lane-Emden equation with index n 3 / 2 = D/2 and 
condition at the origin 9(0) = 1. 

The structure and the stability of polytropic spheres in 
various dimensions of space has been studied by Chavanis 




FIG. 1: Relation between the mass (ordinate) and the cen- 
tral density (abscissa) of box-confined polytropes with index 
n 3 / 2 = D/2 in an appropriate system of coordinates (see 
23, 25] for details). Complete polytropes correspond to the 
ternrinal point in the series of equilibria. The series becomes 
dynamically unstable with respect to the Euler-Poisson sys- 
tem (saddle point of the energy functional) after the turning 
point of mass r\ which appears for D = 4. Thus, classical 
white dwarf stars are stable for D < 4 and unstable for D > 4. 



For 2 < D < 4, the mass M decreases as the radius R 
increases while for D < 2 and for 4 < D < 2(1 + y/2) it 
increases with the radius (see Fig. [5]). 
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FIG. 2: Mass-radius relation for classical white dwarf stars 
in different dimensions of space. The radius is independent 
on mass in D — 2 and the mass is independent on radius in 
D — 4. The dimensional factors Mq and Rq are defined in 
Sec. EH 
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For D > 2, this study exhibits two important 



D 



D 



n 3 = 



D 



n 5 = 



D 



(30) 



According to this study, a classical white dwarf star is 
self-confined (complete) if n 3 / 2 < n 5 , i.e. D < 2(1 + 
y/2) = 4.8284271.... In that case, the density vanishes at 
a finite radius R identified as the radius of the star. On 
the other hand, it is nonlinearly dynamically stable with 
respect to the Euler-Poisson system if n 3 / 2 < 713, and 
linearly unstable otherwise. Therefore, classical white 
dwarf stars are stable for D < 4 and unstable for D > 4 
(see Fig. [T]and Appendix |B|) . For D < 2, classical white 
dwarf stars are always self-confined and stable. Using 
the results of [23], the mass-radius relation for complete 
polytropes with index 713/2 = D/2 in D dimensions is 
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where £1 is such that 6 D / 2 (£i) = 0. Using Eq. (f2"(J)) . 
we find that the mass-radius relation for classical white 
dwarf stars in D dimensions is 
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For D = 4 the mass is independent on the radius and 
given in terms of fundamental constants by 
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We recall that the value of the gravitational constant G 
depends on the dimension of space so that we cannot 
give an explicit value to this limiting mass. The central 
density is related to the radius by 
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For D = 2, the radius is independent on mass and given 
in terms of fundamental constants by 
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Furthermore, for D = 2, the Lane-Emden equation 
with index n 3 / 2 = 1 can be solved analytically, yielding 
0i = ^o(0 and 6 = "0,1 where a ,i = 2.404826... is 
the first zero of Bessel function Jo. We obtain a density 
profile 

p(r) = poMZir/R), (37) 

where the central density is related to the total mass by 
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Finally, for D = 3, using the mass-radius relation (|33|) . 
we find that the average density is related to the to- 
tal mass by p = 2.162 10 6 (M/M©) 2 g/cm 2 (for fj, = 
2.5). Historically, this result was obtained by Chan- 
drasekhar (1931) [8] who first applied the theory of 
polytropic gas spheres with index n = 3/2 to classi- 
cal white dwarf stars. It improves an earlier result 
p = 3.977 10 6 (M/Mq) 2 g/cm 2 obtained by Stoner (1929) 
[3fJ on the basis of his model of stars with uniform den- 
sity (see Sec. |VIT)) . 



V. THE ULTRA-RELATIVISTIC LIMIT 

In the ultra-relativistic limit x 3> 1, we find that the 
equation of state takes the form 
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Therefore, an ultra-relativistic white dwarf star is equiv- 
alent to a polytrope of index 
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FIG. 3: Relation between the mass (ordinate) and the cen- 
tral density (abscissa) of box-confined polytropes with index 
n' 3 — D in an appropriate system of coordinates (see [231 . l25fl 
for details). Complete polytropes correspond to the termi- 
nal point in the series of equilibria. The series becomes dy- 
namically unstable with respect to the Euler-Poisson system 
(saddle point of the energy functional) after the turning point 
of mass 77 which appears for D = 3. Thus, ultra-relativistic 
white dwarf stars are stable for D < 3 and unstable for D > 3. 
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This also directly results from the Chandrasekhar equa- 
tion. For x ^> 1, it reduces to 
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Using Eq. (|4"0"|) , we find that the mass-radius relation for 
ultra-relativistic white dwarf stars in D dimensions is 
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For 3 < D < (3 + VT7)/2, the mass M increases as the 
radius R increases while for 1 < D < 3 it decreases with 
the radius (see Fig. [4j. 



where we have set 9 — <fi and £ = 77. This is the Lane- 
Emden equation with index n' 3 = D. 

Using the results of [23] for D > 2, we deduce that 
an ultra-relativistic white dwarf star is self-confined if 
n' 3 < n 5 , i-e. D < (3 + Vl7)/2 = 3.5615528.... In addi- 
tion, it is nonlinearly dynamically stable with respect to 
the Euler-Poisson system if n' 3 < n^ and linearly unsta- 
ble otherwise. Therefore, ultra-relativistic white dwarf 
stars are stable for D < 3 and unstable for D > 3 (see 
Fig. and Appendix [B]) . For D < 2, ultra-relativistic 
white dwarf stars are self-confined and stable. On the 
other hand, using the results of [23|], the mass-radius re- 
lation for complete polytropes with index n' 3 = D in D 
dimensions is 
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FIG. 4: Mass-radius relation for ultra-relativistic white dwarf 
stars in different dimensions of space. The radius is indepen- 
dent on mass in D — 1 and the mass is independent on radius 
in D — 3 (Chandrasekhar's mass). 



For D = 1, the radius is independent on mass and 
given in terms of fundamental constants by 
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Furthermore, for D = 1, the Lane-Emden equation 
with index n^ = 1 can be solved analytically, yielding 
6*i = cos(£) and £i = 7r/2. We obtain a density profile 
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where the central density is related to the total mass by 
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For D = 3, the mass is independent on radius and given 
in terms of fundamental constants by 
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This is the Chandrasekhar mass 
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Coming back to Eq. (TT9"]) . we can show that for this 
limiting value, the radius R of the configuration tends 
to zero (see Sec. I VI I) . Historically, the existence of 
a maximum mass for relativistic white dwarf stars was 
first published by Anderson (1929) [18] who considered a 
relativistic extension of the model of Stoner (1929) [30| 
for classical white dwarf stars. He obtained a limiting 
mass M = 1.37 10 33 g (for /i = 2.5). The relativis- 
tic treatment of Anderson was criticized and corrected 
by Stoner (1930) [19| who obtained a value of the limit- 
ing mass M = 2.19 10 33 g. The uniform density model 
of Stoner was in turn criticized and corrected by Chan- 
drasekhar (1931) [l(| who applied the theory of poly- 
tropic gas spheres with index n — 3 to relativistic white 
dwarf stars and obtained the value (|5i"j) of the limiting 
mass M — 1.822 10 33 g. It seems that these historical 
details are not well-known because the references to the 
works of Anderson and Stoner progressively disappeared 
from the literature. 



VI. THE GENERAL CASE 



Collecting together the results of Sec. IIII( the mass- 
radius relation for relativistic white dwarf stars in D di- 
mensions can be written in the general case under the 
parametric form 



M 



D(3-D)/2 

Vo 



n(»0, 



R_ 
Rq 



1 



(D-l)/2 



m, (52) 



8A 2 \ D/2 1 
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The mass Mo and the radius Rq can be expressed in terms 
of fundamental constants as 
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We can now obtain the mass-radius curve M — R by the 
following procedure. We fix a value of the parameter yo 
and solve the differential equation (fl9|) with initial con- 
dition ([20]) until the point r\ = rji, determined by Eq. 
(|2"Tj) . at which the density vanishes. The radius and the 
mass of the corresponding configuration are then given 
by Eq. (|52[) . By varying j/Oj we can obtain the full curve 
R(yo) ~ M(yo) parameterized by the value of the central 
density p given by Eq. (|2"3")l . To solve the differential 
equation (|19[) . we need the behavior of <f> at the origin. 
Expanding (f)(rj) in Taylor series and substituting this ex- 
pansion in Eq. (|19[) we obtain for 77 — ► 0: 
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We note in particular that <fi"(0) — —q/D. 

Finally, we give the asymptotic expressions of the 
mass-radius relation. In the classical limit, using Eq. 
(l33l. we obtain 



M 

Mo 



(D-2)/D 



R 

Ro 



4-D 



_ 1 (D-2)/D 

- 2^/2 • 



(59) 



In the ultra-relativistic limit, using Eq. (|4^)l . we get 
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In Figs. [51TT21 we plot the mass-radius relation of rela- 
tivistic white dwarf stars (full line) for different dimen- 
sions of space. The asymptotic relations ([5T)| and (|6"0|) 
valid in the classical (C) and ultra-relativistic (R) lim- 
its are also shown for comparison (dashed line) together 
with the analytical approximation (dotted line) derived 
in Sec. ED 
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FIG. 5: Mass-radius relation in D = 1. The radius increases 
with the mass. The configurations are always stable and there 
exists a maximum radius achieved in the ultra-relativistic 
limit for M —* +oo. 
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FIG. 6: Mass-radius relation in D — 1.5. The configurations 
are always stable and there exists a maximum radius for par- 
tially relativistic distributions. 



VII. ANALYTICAL APPROXIMATION OF THE 
MASS-RADIUS RELATION 



We present here an analytical approximation of the 
mass-radius relation in various dimensions of space based 
on the treatment by Nauenberg [21( in D — 3. This treat- 
ment amounts to considering that the density is uniform 
in the star and the mass-radius relation is obtained by 
minimizing the energy functional (see Appendix [C]) with 
respect to the radius (or to the density) at fixed mass. 
As mentioned in the Introduction, this is similar to the 
simplified model of relativistic white dwarf stars made by 
Stoner [ljjj before Chandrasekhar's treatment [10]. Fol- 
lowing Nauenberg [2l| . we approximate the kinetic en- 




FIG. 7: Mass-radius relation in D — 2. The radius decreases 
with the mass. The configurations are always stable and there 
exists a maximum radius achieved in the classical limit for 
M-> 0. 
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FIG. 8: Mass-radius relation in D — 2.5. The radius decreases 
with the mass. The configurations are always stable. 



ergy K by the form 

K = Nmc 2 { | I 



2 \ 1/2 



(61) 



where N is the number of electrons and p is an average 
over the star of the momentum of the electrons [56fl . We 
assume that it is determined by an appropriate average 
value of the density by the relation 



2<S , d „ D 



(62) 



based on the Pauli exclusion principle. Now, for the den- 
sity, we write 



p = C 



M 



SdRD' 



(63) 
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FIG. 9: Mass-radius relation in D — 3. The radius decreases 
with the mass. The configurations are always stable and there 
exists a limiting mass (Chandrasekhar's mass) achieved in 
the ultra-relativistic limit for R — 0. For M > Mchandra, 
quantum mechanics cannot arrest gravitational collapse. 




FIG. 10: Mass-radius relation in D = 3.3. There exists a 
limiting mass for partially relativistic distributions. Classical 
configurations are stable and the series of relativistic equilib- 
ria becomes unstable after the turning point of mass. Thus, 
highly relativistic configurations in D > 3 cannot be in hy- 
drostatic equilibrium. 



FIG. 11: Mass-radius relation in D = 3.8. For D > 
(3 + v / 17)/2 the ultra-relativistic configurations (equivalent 
to polytropes of index n' 3 — D) are not self-confined anymore 
and a spiral develops in the mass-radius relation. This is 
somehow similar to the classical spiral occurring in the (E, f3) 
plane in the thermodynamics of isothermal self-gravitating 
systems [3U , I32J [331 . |34| and to the spiral occurring in the 
(M, R) plane in the general relativistic treatment of neutron 
stars [351 l3q |. The series of equilibria becomes unstable at 
the first turning point of mass and new modes of instability 
occur at the secondary turning points [37| . 
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where £ is a dimensionless parameter. We also write the 
potential energy in the form 



W = - 



GM 2 



D-2R D ~ 2 



(64) 



where v is another dimensionless parameter. By writing 
Eq. ([64]) . we have assumed that D^2 but we shall see 
that the following results pass to the limit for D — ► 2. 
We introduce two dimensionless variables n and r and 
two fixed constants M* = N^fxH and i?» such that 



M = nlVL 



R = rR* 



(65) 



FIG. 12: Mass-radius relation in D — 4. There exists a limit- 
ing mass achieved in the classical limit for R — » +oo. In fact, 
for D > 4, all the configurations (classical and relativistic) are 
unstable. Quantum mechanics cannot balance gravitational 
attraction even in the classical limit. 



We determine M* and i?* by the requirement that the 
relativity parameter x = p/mc have the form 



AID 



(66) 



and that the potential energy can be written 

,2 



w 



1 9»r U 

-mc ■ N, 



D-2 



D-2 ' 



(67) 



This yields 



R* 



1/2 



h D/2 



DC 

2vS%) G 1 /2 TO (I3-l)/2 c (I3-2)/2 /ijff ' 



(68) 



M, 



1 / D Q\( D - 2 )/ 2 h D(D-2)/2 c D(i-D)/2 1 



v D / 2 \2S 2 D 



m D(D-3)/2 G D/2 (^fl-1 
(69) 

Comparing with Eqs. (|55[) and ([56]) . we find that 



KL 



^ D -^M , **=(£) 



1/2 



Ro- (70) 



Now, the energy E = K + W of the star can be written 

„-D-2 



£ = A^mcrV (1 + x 2 ) 1/2 -l-n 2/ 



D 



D-2 



In the classical limit i < 1, we get 



Edass = JV.mc 



,1+2/B 



2r 2 (L> - 2)r-°- 2 

and in the ultra- relativistic limit i>1. we obtain 



Erelat = N*mC 



,1+l/D 



r (-D - 2) 



,D-2 



(71) 



(72) 



(73) 



We shall consider the energy as a function of the radius 
R, with the mass M fixed. Thus, the mass-radius rela- 
tion will be obtained by minimizing the energy E versus 
x with fixed n. Writing dE/dx — 0, we obtain the equa- 
tions 



x D(i-D)/2 
(l + X 2 ) D / 4 ' 



A/D 



(74) 



(75) 



defining the mass-radius relation in parametric form in 
the framework of the simplified model. Note that x D 
represents the value of the density in units of M*/R®, i.e. 
p = x D M !f /R^ ) . Therefore, Eq. (|T4"|) can be viewed as the 
relation between the mass and the density. Eliminating 
the relativity parameter x (or density) between Eqs. (I74|) 
and (|75|). we explicitly obtain 



VT 



1-r 4 



(D = l) 



(D = 2) 



(76) 



(77) 



(1 - n 4 / 3 ) 1 / 2 

„l/3 



7 3 / 4 



VT 



(D = 3) 



(D = A) 



(78) 



(79) 



In D = 1, there exists a maximum radius r = 1 achieved 
for n — > +oo, and we have the scaling 



1 



;(l-r) 



-1/2 



In D = 2, there exists a maximum radius r 
for n — > 0, and we have the scaling 

n — 4(1 - r). 



(80) 
1 achieved 

(81) 



In D = 3, there exists a maximum mass n = 1 (Chan- 
drasekhar's mass) achieved for r — > 0, and we have the 
scaling 



>-">" 2 - 



(82) 



In D = 4, there exists a maximum mass n = 1 achieved 
for r — > +oo, and we have the scaling 



r ~ (1 — n) 



-1/2 



(83) 



In order to compare Eqs. (|76"|) - (|T9"|) with the exact results, 
we need to estimate the values of the constants C and v. 
This can be done by considering the limiting forms of 
Eqs. ([74]) and ([75]) . In the classical limit i « 1, we 
obtain 



n (D-2)/D r A-D = 1; 

and in the ultra- relativistic limit x 3> 1, we get 



n 



(D-l)/D3-D 



Comparing with Eqs. 



D uj d ~ 2 

, ,D-1- 



and (foT)]) . we find that 



!D-2)/D 
1 W D/2 



2 , 2(D-1)/D' 



(84) 



(85) 



(86) 



where the quantities cod and wd/2 can be deduced from 
the numerical study of the Lane-Emden equation. The 
analytical approximations of the mass-radius relation 
dZ6l)-(IZ9l) are plotted in dotted lines in Figs. 01 
and they give a fair agreement with the exact results 
(full line). Of course, they cannot reproduce the spiral in 
D = 4 which requires the resolution of the full differential 
equation (fi"9]) . 

Let us address the stability of the configurations within 
this simplified analytical model. If we view a white dwarf 
star as a gas of electrons at statistical equilibrium at tem- 
perature T, stable configurations are those that minimize 
the free energy F = E — TS at fixed mass, where E is 
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the total energy (kinetic + potential) and S is the Fermi- 
Dirac entropy (see Appendix[U]). For a completely degen- 
erate gas at T = 0, stable configurations are those that 
minimize the energy E at fixed mass. As we have seen 
previously, the first order variations E'(x) — determine 
the mass-radius relation. Then, the configurations is sta- 
ble if the energy is minimum, i.e. E"(x) > 0. Now, we 
have 



E" 



N+mxrn 

(1+X 2 ) 3 / 2 



[4 - D - (D - 3)x 2 }. 



(87) 



For D < 3, all the solutions are stable while for D > 4, all 
the solutions are unstable. For 3 < D < 4, the solutions 
with x < x c are stable and the solutions with x > x c are 
unstable where 



4-D 
D-3 



1/2 



(88) 



We now show that the onset of instability precisely cor- 
responds to the turning point of mass, i.e. to the point 
where the mass is maximum in the series of equilibria 
M(R). In terms of reduced variables, this corresponds 
to dn/dr = 0, or equivalently n'(x) = 0. Taking the 
logarithmic derivative of Eq. ([71)1 , we have 

dn D , , _. dx D x , ,„„. 

— = — U-D) — -dx, (89) 

n 2 v ' x 2 1 + x 2 ' v ' 

so that the condition n'(x) — yields x = x c . Thus, in- 
stability sets in precisely at the maximum mass as could 
have been directly inferred from the turning point crite- 
rion |37j . 

It is also useful to discuss the classical and the ultra- 
relativistic limits specifically. In the classical limit, using 
Eq. flTJ), we get 



F" 

class 



N*mc 2 n(A- D), 



(90) 



so that a classical white dwarf star is stable for D < 4 
and unstable for D > 4. In the ultra-relativistic limit, 
using Eq. (f73|) . we get 



E. 



N*mc n 



dat 



(3-D), 



(91) 



so that an ultra-relativistic white dwarf star is stable for 
D < 3 and unstable for D > 3. This returns the results 
obtained in Sees. IIVI and IVl 

Finally, considering the function E(R) given by Eq. 
(I71[) and assuming that the system evolves so as to mini- 
mize its energy (this requires some source of dissipation), 
we have the following results [53]: for D < 3, there ex- 
ists a stable equilibrium state (global minimum of E) 
with radius R for all mass M . For D = 3, there exists 
a maximum mass M c so that: (i) for M < M c , there 
exists a stable equilibrium state (global minimum of E) 
with radius R > (ii) for M = M c , the system collapses 
to a point (R = 0) but its energy remains finite (lower 



bound) (iii) for M > M c , the system collapses to a point 
R — » and E — > — oo. For 3 < D < 4, there exists a 
maximum mass M c so that: (i) for M < M c , there ex- 
ists a metastable equilibrium state (local minimum of E) 
with x < x c and an unstable equilibrium state (global 
maximum of E) with x > x c . The system can either 
reach the metastable state or collapse to a point (R —* 0, 
E — > — oo); the choice probably depends on a notion of 
basin of attraction (ii) for M > M c , the system collapses 
to a point (R — * 0, E — > — oo). For D — 4, there exists 
a critical mass M c so that: (i) for M < M c , there exists 
an unstable equilibrium state (global maximum of E) so 
the system either collapses (R — ► and E — ► — oo) or 
evaporates (R — > +oo and E — > 0) (ii) for M > M c , the 
system collapses to a point (R — > and E — > — oo). For 
D > 4, there exists an unstable equilibrium state (global 
maximum of i?) for all mass M so the system either col- 
lapses (R — + and E — ► — oo) or evaporates (i? — ► +oo 
and£^0). 



VIII. CONCLUSION 

"Why is our universe three dimensional? Does the di- 
mension D = 3 pZay a special role among other space 
dimensions?" 

Several scientists have examined the role played by the 
dimension of space in determining the form of the laws 
of physics. This question goes back to Ptolemy who ar- 
gues in his treatise On dimensionality that no more than 
three spatial dimensions are possible in Nature. In the 
18 th century, Kant realizes the deep connection between 
the inverse square law of gravitation and the existence 
of three spatial dimensions. Interestingly, he regards the 
three spatial dimensions as a consequence of Newton's 
inverse square law rather than the converse. In the twen- 
tieth century, Ehrenfcst [26], in a paper called "In what 
way does it become manifest in the fundamental laws of 
physics that space has three dimensions?" argues that 
planetary orbits, atoms and molecules would be unstable 
in a space of dimension D > 4. This idea has been fol- 
lowed more recently by Gurevich & Mostepanenko [38| 
who argue that if the universe is made of metagalaxies 
with various number of dimensions, atomic matter and 
life are possible only in 3-dimensional space. Other in- 
vestigations on dimensionality are reviewed in the paper 
of Barrow [39J . We have found that the relativistic self- 
gravitating Fermi gas at T — (a white dwarf star) pos- 
sesses a rich structure as a function of the dimension of 
space. We have exhibited several characteristic dimen- 
sions D = 1, D = 2, D = 3, D = (3 + VT7)/2, D = 4 
and D = 2(1 + V2). For D < 3, there exists stable con- 
figurations for any value of the mass. For D = 3, the se- 
quence of equilibrium configurations is stable but there 
exists a maximum mass (Chandrasekhar's limit) above 
which there is no equilibrium state. For 3 < D < 4, 
the sequence of equilibrium configurations is stable for 
classical white dwarf stars but it becomes unstable for 
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relativistic white dwarf stars with high density after the 
turning point of mass. Therefore, the dimension D = 3 
is special because it is the largest dimension at which the 
sequence of equilibrium configurations is stable all the 
way long; for D > 3, a turning point of mass appears 
so that ultra-relativistic white dwarf stars become unsta- 
ble. Therefore, the dimension D — 3 is marginal in that 
respect. Finally, for D > 4, the whole sequence of equi- 
librium configurations is unstable. Therefore, as already 
noted in our previous papers [23|, [24], [25| , the dimension 
D = 4 is critical because at that dimension quantum 
mechanics cannot stabilize matter against gravitational 
collapse, even in the classical regime, contrary to the sit- 
uation in D = 3 [J, [jj Q. Interestingly, this result is 
similar to that of Ehrenfest [2(| although it applies to 
white dwarf stars instead of atoms. 

Our exact description of D-dimensional white dwarf 
stars based on Chandrasekhar's seminal paper [l3j shows 
that relativistic white dwarf stars become unstable in 
D > 3 and that classical white dwarf stars become un- 
stable in D > 4. Therefore, for D > 4, a self-gravitating 
Fermi gas forms a black hole or evaporates (see Appendix 
|B"|) . These conclusions have also been reached by Bech- 
hoefer & Chabricr 28] on the basis of simple dimensional 
analysis. This suggests that a D > 4 universe is not vi- 
able (see also Appendix [DJ and gives insight why our 
universe in apparently three dimensional. We note that 
extra-dimensions can appear at the micro-scale, an idea 
originating from Kaluza-Klein theory. This idea took a 
renaissance in modern theories of grand unification which 
are formulated in higher-dimensional spaces [58[ • Our ap- 
proach shows that already at a simple level, the coupling 
between Newton's equations (gravitation), Fermi-Dirac 
statistics (quantum mechanics) and special relativity re- 
veals a rich structure as a function of D. In this respect, 
it is interesting to note that the critical masses (|3~4"|) (|5D|) 
and radii (|3l)|) Wf\ that we have found occur for sim- 
ple integer dimensions D = 1, 2, 3 and 4, which was not 
granted a priori. 

It is interesting to develop a parallel between the mass- 
radius relation M(R) of white dwarf stars and the caloric 
curve T(E) giving the temperature as a function of the 
energy in the thermodynamics of self-gravitating systems 
[3ll 132 . 1331 [34[ . In this analogy, the Chandrasekhar mass 
in D = 3 is the counterpart of the critical temperature 
for isothermal systems in D = 2 (in both cases, the equi- 
librium density profile is a Dirac peak containing all the 
mass) 59]. On the other hand, for D > 3, the mass- 
radius relation for white dwarf stars exhibits turning 
points, and even a spiraling behavior for D > i(3 + VT7), 
which is similar to the spiraling behavior of the caloric 
curve for isothermal systems in D = 3. In this analogy, 
the maximum mass, corresponding to a critical value of 
the central density which parameterizes the series of equi- 
libria, is the counterpart of the Antonov energy (in the 
microcanonical ensemble) [3l| or of the Emden tempera- 
ture (in the canonical ensemble) [34( . The series of equi- 
libria becomes unstable after this turning point. In addi- 



tion, there is no equilibrium state above this maximum 
mass, or below the minimum energy or minimum tem- 
perature in the thermodynamical problem. In that case, 
the system is expected to undergo gravitational collapse. 

As a last comment (notified by the referee), it should 
be emphasized that the conclusions reached in this pa- 
per concerning dimensionality implicitly assume that the 
laws of physics that we know remain the same in a uni- 
verse of arbitrary dimension D. This is of course not 
granted at all. There may be a new cosmological theory, 
a new theory of star formation and stellar evolution in 
higher dimensions. We also emphasize that our approach 
does not take into account general relativistic effects that 
can sensibly modify the results [43 ]. 
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APPENDIX A: THE EQUATION OF STATE 

We provide here the explicit expression of the function 
fix) defined by Eq. (TTUJ) for different dimensions of space 
D = 1, 2, 3 and 4, respectively: 



f{x) = 4a;(l + x 2 ) 1/2 - 4 sink -1 5 



(Al) 



/(*) = ^ + ^ 2 -2)(l + * 2 ) 1/2 , (A2) 



f(x) = x(2x 2 - 3)(1 + x 2 ) 1 ' 2 + 3 sinh -1 ^ (A3) 



/(*) = ~ + ^(3x 4 - Ax 2 + 8)(1 + x 2 fl 2 . (A4) 
15 15 



APPENDIX B: STABILITY CRITERIA FOR 
POLYTROPIC SPHERES IN D DIMENSIONS 

We generalize in D dimensions the usual stability cri- 
teria for polytropic gaseous spheres, and apply them to 
classical and ultra-relativistic white dwarf stars. 



1. The Euler-Poisson system 



Let us consider the Euler-Poisson system [43] : 

^+V-( /O u) = 0, (Bl) 



o 1 

— + (u-V)u= — VP-V$, (B2) 

at p 



A$ = S D Gp, (B3) 
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describing the dynamical evolution of a barotropic gas 
with an equation of state P = P(p). The Euler-Poisson 
system conserves the mass M and the energy [251 ]: 



W 



j P f '^T d P' dr +l J Pfdv+ l - J pu 2 dr. (B4) 



In the following, we shall essentially consider the case of 
a polytropic equation of state P(p) = Kp 1 where 7 = 
1 + 1/n. For 7^1, the energy functional W is given by 

W = — j- /Vdr+i [p$dr+^ f pu 2 dr. (B5) 
It can be rewritten 
W-- 



7 



— / Pdr + - / p<5>dr +- pu 2 dr. (B6) 



For an isothermal equation of state P = pkgT '/to, corre- 
sponding to 7 — ► 1 or n — ► +00, we have 



W = k B T 



— In — dr ' 



1 



., p<Mr+- / pu 2 dr. (B7) 



For a general polytropic equation of state P = Kp 1 with 
arbitrary 7, it is convenient to write the energy of the 
gas in the form 



W 



K 



7 



-J{p^-p)dv+ l -Jp^dv+ l -Jpn 2 dv. (B8) 



We have added a constant term — -^j J pdr proportional 
to the total mass (which is a conserved quantity) so as to 
recover the energy (|F37|) of an isothermal gas in the limit 

7->l M- 



2. The eigenvalue equation 

We first consider the linear dynamical stability of 
a polytropic star described by the Euler-Poisson sys- 
tem and generalize the app roach developed in Chavanis 
[H [H, M, [H, |H |H, 53 ■ We consider a steady solu- 
tion of the Euler-Poisson system satisfying u = and 
the condition of hydrostatic balance VP + pV^ = 0. For 
a polytropic equation of state P = Kp 1+1 / n , the equi- 
librium density profile is solution of the D-dimensional 
Lane-Emden equation [231 ] - We consider complete poly- 
tropes such that the density vanishes at a finite radius R. 
For n 7^ n 3 , there exists a unique steady state for any 
mass M. The mass-radius relation is [231 ]: 

M {n-l)/n R (D-2)(n 3 -n)/n _ K ( 1 + n ) ,(n-l)/n mgX 



GS* 



u 



where 



. e (n + l)/ (n -l)^ (ei) 



(BIO) 



is a constant given in terms of the solution 9 n (£) of the 
Lane-Emden equation of index n in D dimensions. For 



the critical index n = n^, steady state solutions exist 
only for a unique value of the mass 



M„ 



K(l 



GS 



l/n 
D 



n/(n-l) 



UJ,, 



(Bll) 



These solutions have the same mass M c but an arbitrary 
radius R. We already anticipate that the index n = n.3 
will correspond to a case of marginal stability separating 
stable and unstable solutions. 

Linearizing the equations of motion (|B1|) - (|B3[) around 
a stationary solution in hydrostatic balance and writing 
the perturbation in the form 5p ~ e xt , we obtain after 
some calculations [3J, H|| H3] the eigenvalue equation 



P'{p) dq 
dr ySopr ^ 1 dr 



Gq_ 
„D-1 



X 2 



S D pr ] 



j^Tl, (B12) 



where we have restricted ourselves to spherically sym- 
metric perturbations and defined q(r) = J 5pSDr D ~ 1 dr. 
For a polytropic gas with an equation of state P = Kp" f , 
the foregoing equation becomes [45| | : 



K~f- 



P 



,7-2 



dq 



dr \ Snr D 1 dr 



Gg 

r D-l 



X 2 



S D pr D - 



-q. (B13) 



The polytrope is stable if A 2 < (yielding oscillatory 
modes with pulsation u — \J— A 2 ) and unstable if A 2 > 
(yielding exponentially growing modes with growth rate 
7 = vA 2 ). Considering the point of marginal stability 
(A = 0) and introducing the Emden variables [231 [29J. l4a| . 
Eq. (|5l3|) reduces to 



d_ 
dr 



1 -™ dF 



nF 



$ 



D-l 



= 0. 



This equation has the exact solution [23|, |45( : 



F(0 = ci 



n — 1 



(B14) 



(B15) 



The point of marginal stability is then determined by 
the boundary conditions. One can show [45|, [47| that the 
velocity perturbation is given by Su = —Xq/(SDpr D ^ 1 )- 
Therefore, if the density of the configuration vanishes at 
r = R, one must have q(R) = to avoid unphysical diver- 
gences. Thus, if £1 denotes the value of the normalized 
radius R of the star [H, [H |Zi| such that 6>(£i) = 0, the 
natural boundary condition for the eigenvalue equation 
(IB13P is F(£i) = 0. Substituting this condition in Eq. 
(IB15[) . we obtain the critical index 



D 



D-2 



n 3 , 



(B16) 



corresponding to a marginally stable gaseous polytrope 
(A = 0). Our method gives the form (|B15j) of the neu- 
tral perturbation Sp at the critical index n = n 3 . We 
conclude that the infinite family of steady state solutions 
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with equal mass M c and arbitrary radius R that exists 
for n — n 3 is marginally stable. 

We now consider the nonlinear dynamical stability 
problem. Since the Euler-Poisson system (|B1[) - (|B3[) con- 
serves the mass M and the energy W, a maximum or 
a minimum of the energy functional W[p, u] at fixed 
mass M[p] — M determines a steady state of the Euler- 
Poisson system that is formally nonlinearly dynamically 
stable [25] . Because of the presence of the kinetic term 
0[u] = (1/2) J pu 2 dr, the functional W[p, u] has no ab- 
solute maximum. Thus, we need to investigate the pos- 
sible existence of a minimum of W[p, u] at fixed mass 
M[p] = M. A barotropic star that minimizes the energy 
functional W at fixed mass M is nonlinearly dynamically 
stable with respect to the Euler-Poisson system. The 
cancellation of the first order variations SW — aSM = 0, 
where a is a Lagrange multiplier, yields u = and the 
condition of hydrostatic balance VP + pV$ = 0. Then, 
the condition of nonlinear dynamical stability is 



5 2 W 



P'(P) 
2p 



(Spfdr 



5p8$dv > 0, (B17) 



for all perturbations that conserve mass, i.e. J Spdr — 
0. After some calculations [25j, this can be put in the 
quadratic form 



<rw 



d 
dr 



P'{p) dq 



S D pr 



D-l 



dr 



Gq_ 

r.D-1 



q dr. 
(B18) 



We are led therefore to consider the eigenvalue problem 



d_ 
dr 



P'{p) 



S D pr 



D-l 



dr 



G 



n D-l 



qx(r) = Xqx(r). (B19) 



If all the eigenvalues A are negative, then 8 2 W > and 
the configuration is a minimum of W at fixed mass. This 
implies that it is nonlinearly dynamically stable. If at 
least one eigenvalue A is positive, the configuration is a 
saddle point of W and the star is dynamically unstable. 
The marginal case is when the largest eigenvalue A is 
equal to zero. Now, for A = 0, Eqs. ljB19|) and (|B12|) 
coincide. This implies that the conditions of linear and 
nonlinear dynamical stability are the same. In the case 
of polytropic stars, the case of marginal stability corre- 
sponds to the critical index n = n 3 . At that index, the 
equilibrium configurations with mass M c and radius R 
all have the same value of energy W = (sec Eq. ([B43P 
later). This is therefore a very degenerate situation. 

The nonlinear dynamical stability of gaseous poly- 
tropes can also be investigated by plotting the series 
of equilibria M(po) (mass vs central density) of box- 
confined configurations and using the turning point ar- 
gument of Poincare [23j, [25J, |45j . We can thus determine 
whether the last point on the series of equilibria, which 
corresponds to a complete polytrope whose density van- 
ishes precisely at the box radius, is stable or not. For 
D < 2, there is no turning point of mass (we restrict our- 
selves to n > 0), implying that the gaseous polytropes 



are always stable. For D > 2, a turning point of mass 
M(po) appears precisely for n — n 3 (see Fig. 5 of [23]]). 
This method shows that, for D > 2, complete polytropes 
with n < ri3 are nonlinearly dynamically stable (they are 
minima of W at fixed mass M) while complete polytropes 
with n > 7J3 are dynamically unstable (they are saddle 
points of W at fixed mass M). Complete polytropes with 
n — 713 and M = M c are marginally stable. 



3. The Ledoux criterion 

We can also investigate the linear dynamical stability 
of gaseous polytropic spheres by using the method intro- 
duced by Eddington [48] and Ledoux [49| . If we introduce 
the radial displacement 



£(r) 



6 a 



Sdpt 1 



(B20) 



we can rewrite the eigenvalue equation (|B13|) in the form 
[H,ll: 



d 

dr 



p 7r D+l<% \ , ,.D 



dr 



dP 



+ r u {D 1 + 2-2D)^—i 
dr 



X 2 pr 



D+l 



(B21) 



This is the Eddington equation of pulsations, which has 
been written here in the form of a Sturm-Liouville prob- 
lem. It must be supplemented by the boundary condi- 
tions 



Sr = £r = 0, in r = 0, 



dP = A 7 P [ D^ + r$ 
dr 



0. 



(B22) 



R, (B23) 



where dP/dt — dSP/dt + SudP/dr is the Lagrangian 
derivative of the pressure. Since P = at the surface 
of the star, it is sufficient to demand that £ and d^/dr 
be finite in r = R. Multiplying Eq. (|B21[) by £ and 
integrating between and R, we obtain 



pr D+i edr 



R dr\p V ^(f' 



t2„D 



£ 2 r u {Dj + 2-2D) 



dr 

dP 
dr 



(B24) 



The system is linearly dynamically stable if A 2 < and 
unstable otherwise. Since dP/dr < 0, a sufficient condi- 
tion of stability is D7 + 2 - 2D > 0, i.e. 



7 > 74/3 



2(J>-1) 
D 



1 1 

n ri3 



D 



D 



(B25) 



It can be shown furthermore (see below) that the system 
is unstable if 7 < 74/3 so that the criterion (|B25|) is a 
necessary and sufficient condition of dynamical stability 
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(the case 7 = 74/3 is marginal). In terms of the index n, 
a complete polytrope is stable with respect to the Eulcr- 
Poisson system in D = 1 for n > and for n < — 1, in 
D = 2 forn > 0, in D = 3 for < n < 3 and in L> > 2 
for < n < n^. 

From the theory of Sturm-Liouville problems, it is 
known that expression (|B24[) , which can be written A 2 = 
/[£], forms the basis of a variational principle. The func- 
tion £(r) which maximizes the functional /[£] is the fun- 
damental eigenfunction and the maximum value of this 
functional gives the fundamental eigenvalue A 2 . Further- 
more, any trial function under- estimates the value of A 2 
so this variational principle may prove the existence of in- 
stability but can only give approximate information con- 
cerning stability. As shown by Ledoux & Pekeris [49[ , we 
can get a good approximation of the fundamental eigen- 
value by taking £(r) to be a constant (note that £ = Cst., 
i.e. Sr oc r, is the exact solution of the Sturm-Liouville 
equation (|B21|) at the point of marginal stability A = 
for a polytropic equation of state) . For the trial function 
£ = Cst., expression (|B24|) gives 



A^ 



pr 2 dr = (D7 + 2 -2D) 



dP , 
r——dr. 

dr 



(B26) 



According to the condition of hydrostatic balance, 

dP d<f> , 

~T = -P~r, B27 

dr dr 



we have 



dP 

r——dr = — I pr ■ V<&dr, 

dr 



(B28) 



is the potential energy and 



$ 



G 



P(r') 



D-2 I r-r 



r'\D-2 



dr', 



(B35) 



is the gravitational potential. Therefore, we can rewrite 
Eq. (|Bl6|) in the form 



where 



A 2 = (.D7 + 2 - 2D) 



1=1 pr 2 dr, 



I ' 



(B36) 



(B37) 



is the moment of inertia. Since Wu < 0, we conclude 
that the system is unstable if D7 + 2 — 2D < 0, which 
completes the proof above. On the other hand, Eq. (|B36|) 
provides an estimate of the pulsation period u> = y— A 2 
when the system is stable. This is the D-dimensional 
generalization of the Ledoux stability criterion. 

Note finally that, for a spherically symmetric system, 
using the Gauss theorem 



^ GM(r) 



r D-l 



(B38) 



we get 
W U = -S D G I p(r)M(r)rdr 



GM(r) IW/ . 



(B39) 



where we recognize the Virial 



W,, 



pr ■ V<I>dr. 



(B29) 



Inserting 



This expression may be useful to calculate the potential 
energy. 



4. Virial theorem and Poincare argument 



™ = G I P(r')T^Tndr\ 



(B30) 



in Eq. (|B29j) . interchanging the dummy variables r and 
r' and taking the half-sum of the resulting expressions, 
we get 

W ii = -^f , Pir)p( il drdr'. (B31) 



2 J \r-r'\ D - 2 
Therefore, the Virial can be written 

Wu = (D-2)W, (D£2), (B32) 



The Virial theorem associated with the barotropic 
Euler-Poisson system pl jl -pg ]! is (H,^: 



Id 2 / 
2 dt 2 



= 29 + n + Wu 



(B40) 



where 6 = (1/2) J pu 2 dr is the kinetic energy of the 
macroscopic motion and H = D J Pdr. For a polytropic 
equation of state, the energy functional (|B4[) . which is a 
conserved quantity, can be written 



w = u + w + e, 



(B41) 



where 



Wu = -^~, (D = 2), 



W=- J p$dr, 



(B33) 



(B34) 



where 



U = 



7-1 



Pdr, 



(B42) 



is the internal energy. We note that H = D(j — 1)U = 
(D/n)U for a polytrope. At equilibrium, I = = 0, we 
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get (D/n)U + W u = and W = U + W. For D = 2, this 
implies that U = nGM 2 /A. For D^2, this implies that 



W=[1--)W, 



(B43) 



For W > 0, we find that /(£) -4 +oo for t — > +oo so 
that the system evaporates. Alternatively, for W < 0, 
we find that /(£) — > in a finite time, so that the system 
collapses and forms a Dirac peak in a finite time. If we 
consider the dimension D = 3, Eq. (|B49[) reduces to 



where we recall that 



W 



G 



2(D-2) 



p(r)p(r>) 



r — r 



./ID— 2 



drdr' . 



(B44) 



For D > 2 and n > 0, the system is dynamically stable 
if W < and unstable otherwise (Poincare argument 
|li|). Since W < 0, we find from Eq. IJB43|l that the 
polytropic star is stable if and only if n < 713. At the 
point of marginal stability n = 713, where several steady 
configurations exist with the same mass M c given by Eq. 
(jBlip and an arbitrary radius R, the energy of these 
configurations is W = for all R. 

Using W = (n/D)H + W + Q and eliminating the pres- 
sure term in Eq. (|B40[) , the Virial theorem for D/2 can 
be put in the form 



2 dt 2 V n 



D 

e + — w- 

n 



D - 2 ) W. (B45) 

71 



Alternatively, eliminating the kinetic energy in Eq. 
(|B40p . the Virial theorem for D ^ 2 can be written 



2~dP 



2W 



D 

2 

n 



U + (D- A)W. (B46) 



For 11 = n.3/2, corresponding to classical white dwarf 
stars, the Virial theorem becomes 



2dt 



^=™ + {D-A)W. 



(B47) 



For n — rig, corresponding to relativistic white dwarf 
stars, we have 



l£l_ 

2~di? 



= Q + W+(D- 3)W. 



(B48) 



Finally, considering the polytropic index n = n^ of 
marginal stability, the Virial theorem (|B45[) takes the 
form 



1 - d J = ^-D)B + (D-2)W. 

If we consider the dimension D = 4, it reduces to 
d 2 I 



dt 2 



4VV. 



(B49) 



(B50) 



This equation describes the case of classical white dwarf 
stars at the critical dimension D = 4 where 113/2 = TI3 = 
2. It can be integrated into 



I{t) = 2Wt 2 + dt + C 2 . 



(B51) 



kg-*™ 



(B52) 



This equation describes the case of relativistic white 
dwarf stars at the critical dimension D = 3 where 
n 3 = n 3 = 3. For W > 0, we find that I(t) -v +00 
for £ — > +00 indicating that the system evaporates. 



5. Dimensional analysis 

Finally, we show that the instability criterion for poly- 
tropic stars can be obtained from simple dimensional 
analysis. We shall approximate the internal energy (|B42|) 
and the potential energy (|B34j) by 



U = 



W 



K( 



7-1 \R D 



M 



R l 



GM 2 



D-2R D - 2 ' 



(B53) 



(B54) 



where M is the total mass of the configuration and R 
its radius {y and £ are dimensionless parameters). For 
homogeneous spheres, the values of £ and v are given by 
Eq. (|B70|) . With these expressions, the energy functional 
(IB5I) becomes 



W 



7-1 



^) , « D 



GM 2 



D - 2 i?^- 2 



6. (B55) 



We now need to minimize this functional with respect to 
R for a given mass M (a minimum necessarily requires 
O = 0). We first look for the existence of critical points 
(extrema). The cancellation of the first order variations 



dW 
dR 







-KCDM^R^ 1 -^- 1 + GvM 2 R}- D , (B56) 



yields the mass-radius relation 

_/\ / f("- 1 )/"^(D-2)(n 3 -n)/n 



KC,D 
Gv ■ 



(B57) 



This relation determines the radius R of the star as a 
function of its mass M. This expression is consistent with 
the exact mass-radius relation (|B9[) deduced from the 
Lane-Emden equation (of course, our simple approach 
can only model compact density profiles corresponding 
to n < n 5 for D > 2) [23|]. For n ^ 713, there is one, and 
only one, extremum of W(R) for each mass M. In order 
to have a true minimum, we need to impose 

^ = -GvDM 2 R- D (— - -) > 0. (B58) 

dR z \ 71 3 n J 
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Therefore, the system is stable if 



1 1 
- > — : 

n 773 



(B59) 



and unstable otherwise. This simple dimensional analysis 
returns the exact stability criterion (|B25|) . 

We can be a little more precise. For 1/n > 1/773, the 
functional W(R) has a global minimum reached for a 
finite, and non-zero, value of R. This solution is stable. 
For 1/n < 1/713, the functional W(R) has an unstable 
global maximum: for D < 2, W — > — 00 when R — > +00 
(evaporation) and W — v when i? — v (collapse); for 
Z? = 2, W -> -co when i? -»■ and i? -+ +00; for D > 2 
and n < 0, W — > —00 when R — ► and i? — ► +00; 
for 13 > 2 and 71 > 773, W — v —00 when i? — > and 
W — > when i? — > +00. Since we have constructed 
a particular configuration (homogeneous sphere) which 
makes the energy W diverge to —00, the above arguments 
prove that the exact functional W[p, u] given by Eq. (|B4|) 
has no absolute minimum at fixed mass for 1/n < 1/773. 
Since the functional W[p, u] has only one critical point 
(cancelling the first variations) at fixed mass, we conclude 
that, for 1/n < 1/773, this critical point is a saddle point. 
Therefore, there is no stable steady state of polytropic 
spheres for 1/n < 1/773. For the critical index 77 = 773, 
the relation (|B57[) shows the existence of a critical mass 



M r 



( K(D \ D/2 



The functional (|B55[) can be rewritten 



W = K( 



D M 2{ - D -^I D 
D-2 R D - 2 



M 

M~ c 



2/D' 



(B60) 



■0. 



(B61) 

For M = M c , W(R) = for all R (at equilibrium 9 = 0). 
This determines an infinite family of solutions with the 
same mass M c and different radii. These solutions have 
the same energy and are marginally stable. This returns 
the result of the exact model where the configurations are 
solution of the Lane-Emden equation (see Sec. IB 2[) . For 
M < M c , the function W(R) is monotonically decreasing 
with R (so the system tends to evaporate): for 13 < 2, 
W(R) goes from to —00 and for D > 2, W(R) goes 
from +cx) to 0. For M > M c , the function W{R) is 
monotonically increasing with R (so the system tends to 
collapse): for D < 2, W(R) goes from to +00 and for 
D > 2, W{R) goes from -00 to 0. 

We can obtain a simple dynamical model by using the 
Virial theorem (|B40[) . At equilibrium (I = = 0), we 
have 



— U + (D - 2)W = 0. 

77 



(B62) 



Inserting the expressions (IB53|) and (jB54[) in Eq. (IB62|I . 
we recover the mass-radius relation (JB57I) . For 77 ^ 773 
and for any given mass M, there is only one steady state. 



Its radius Ro given by Eq. (|B57[) and its energy Wo is 
given by Eq. (|B43|) . It corresponds to the extremum 
value of W(R). Estimating the moment of inertia by 



aMR 1 



(B63) 



and inserting the expressions (|B53|) and (|B54[) in Eq. 
(|B46|) . we obtain 



laM^-^- = 2W+(D- 2n)KCM 1+1 ' n R- D l'' 
2 dt z 

v(D - 4) GM 2 



D-2 R D ~ 2 ' 



(B64) 



This equation determines the evolution of the radius of 
the star for a fixed mass M and a fixed energy W. The 
evolution of the kinetic energy is then given by = W — 
U—W. The solution of Eq. (|B64[) depends on two control 
parameters M and W and on the initial condition i?(0) 
and i?(0). We shall consider the case where W is equal 
to the value Wo corresponding to the steady state, such 
that the r.h.s. of the above equation is equal to zero at 
equilibrium. Then, we can rewrite Eq. (IB64J) as 



1 n d 2 R 2 
2 aM ^ 



(D - 2n)KCM l+l l n {R- D ' n - R Q D/ ' 

u(D-A) (GM 2 GM 2 
D-2 \R D -- 



) 



, (B65) 
R 



where Rq is the radius of the star at equilibrium. Since 
the dynamics is non-dissipative [60] , the system does not 
evolve towards the minimum of W(R) (unless R — Rq ini- 
tially). The system can either oscillate around the min- 
imum (stable) or evolve away from it (unstable). Con- 
sidering small perturbations around equilibrium, writing 
R = Rq(1 + e) with e -c 1, and linearizing the foregoing 
equation, we obtain 



|l + (D7 + 2 - 2D) !^™ „, 



(it 2 



aRh 



(B66) 



where we have used the mass-radius relation (|B57|) to 
simplify the expression. This is the equation for a har- 
monic oscillator with pulsation 



lu 2 = (D-f + 2- 2D 



vGM 
aR D ' 



(B67) 



The star is stable if to 2 > and unstable otherwise. This 
returns the exact stability criterion (|B25|) . Furthermore, 
using Eqs. (|B54|) and (|B63|) . the pulsation can be rewrit- 
ten in the form 



w 



-{D 1 + 2-2D) 



W„ 



(B68) 



which exactly coincides with the Ledoux formula (|B36[) . 
Therefore, our simple dimensional model allows to ob- 
tain a lot of interesting results. The case of arbitrary 
perturbations around equilibrium will be considered in a 
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future work. Finally, we note that for n = 713, Eq. (JB64I) 
becomes 



1 , r d 2 R 2 
2 aM ^ 



= 2W + K( 



D(D - 4) M 2 ^- 1 )/- 



D-2 
1- 



x 



R D - 2 



(B69) 



and Eq. (|B65|) corresponds to W = and M = M c yield- 
ing d 2 R 2 /dt 2 — 0. Hence, R 2 — C\t + C2 corresponding 
to a marginal evolution. 

In our dimensional analysis, the constants £, v and a 
are dimensionless parameters that could be chosen to fit 
the exact results at best. Alternatively, we can try to ob- 
tain quantitative predictions by calculating their values 
for homogeneous spheres. This yields 



c 



D 



l/n 



D 



v = a = 



D 



(B70) 



Note that the potential energy of a homogeneous sphere 
in D dimensions can be easily calculated from Eq. (1B39J) . 
On the other hand, in D = 2, a direct calculation gives 
W = (1/2)GM 2 \n(R/L)-GM 2 /8 where L is a reference 
radius where $(-L) = 0. All the results given above pass 
to the limit D — ► 2 provided that we take v = 1/2. 
Using Eq. (|B70[) . we find that the approximate value of 
the pulsation (IB67|) becomes 



= (Z>y + 2 - 2D) 



GM 



(B71) 



On the other hand, if we compare the approximate mass- 
radius relation (|B57|) with the exact mass-radius relation 
(IB9|) . we obtain an estimate of the constant u> n in the 
form 



approx. 



D 



1 



n/(rs-l) 



£)l/("-l). 



(B72) 



This has to be compared with the exact value (|B10|) given 
in terms of the solution 9 n (£) of the Emden equation of 
index n in D dimensions [23J ]. Let us consider the case 
D = 3. For n = 3/2, corresponding to classical white 
dwarf stars, we find Wg^™^ = 72 instead of the exact 
value W3/2 = 132.3843.... For n = 3, corresponding to 
relativistic white dwarf stars, we find UJ 'W'P rox - = 2.42... 
instead of the exact value w 3 / 2 = 2.01824... This suggests 
that the homogeneous star model will provide a fair de- 
scription of relativisitic white dwarf stars and a poorer 
description of classical white dwarf stars. We shall come 
back to these different issues (static and dynamics) in a 
future work. 



APPENDIX C: ENERGY FUNCTIONALS 

In this Appendix, we show that the condition of ther- 
modynamical stability in the canonical ensemble is equiv- 
alent to the condition of nonlinear dynamical stability 



with respect to the barotropic Euler-Poisson system. We 
apply this result to white dwarf stars. 



1. Energy of a barotropic gas 

We consider a barotropic gas with an equation of state 
P = P(p) described by the Euler-Poisson system [43| . 
We introduce the energy functional 



w = 



I pf ^r d P' dr + \J P® dr + \J P^ dv - 



(Cl) 

The first term Wi is the work —P(p)d(l/p) done in 
compressing the system from infinite dilution, the sec- 
ond term W is the gravitational energy and the third 
term O is the kinetic energy associated with the mean 
motion. For a gas in local thermodynamic equilibrium, 
the equation of state is P = P(p,T) or P = P(p,s) 
and the first law of thermodynamics can be written 
d(u/p) = —Pd(l/p) + Td(s/p) where s is the entropy 
density and u is the density of internal energy. For a 
gas without interaction (apart from the long-range grav- 
itational attraction), the internal energy is equal to the 
kinetic energy. There are two important cases where the 
gas is barotropic. The first case is when Td(s/p) = 0. 
This concerns either adiabatic (or isentropic) fluids so 
that d(s/p) — or fluids at zero temperature so that 
T = 0. When Td(s/p) = 0, the first law of thermody- 
namics reduces to d(u/p) = —Pd(l/p). It can be inte- 
grated into u = pJ p [P(p')/p' 2 ]dp'. Then, the work Wi 
done by the pressure force (first term in Eq. (|C1[) ) coin- 
cides with the internal energy U of the gas. In that case, 
we get Wi = U and the energy functional (|Clj) can be 
written 



w = u + w + e = E. 



(C2) 



Thus, at T = or for an adiabatic evolution, the 
total energy of the gas E is conserved by the Euler- 
Poisson system (since W is conserved). Alternatively, 
for an isothermal gas dT = 0, the first law of thermo- 
dynamics d(u/p) = —Pd(l/p) + Td(s/ p) can be written 
d((u — Ts)/p) = —Pd(l/p). It can be integrated into 
u-Ts = pf p [P{p')/p' 2 ]dp'. Therefore, the work Wi 
done by the pressure force (first term in Eq. (|C1|) ) co- 
incides with the free energy U — TS of the gas. In that 
case, we get Wi = U — TS and the energy functional 
(IC1D can be written 



W = U-TS + W + B = E-TS = F. 



(C3) 



Thus, for an isothermal evolution, the free energy of the 
system F is conserved by the Euler-Poisson system (since 
W is conserved). At T = 0, we recover the conservation 
of the energy E. 

Let us apply these results to white dwarf stars. We can 
view a white dwarf star as a barotropic gas described by 
an equation of state P = P(p). According to the dis- 
cussion of Appendix IB 21 (see also [25j|). it is nonlinearly 
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dynamically stable with respect to the Euler-Poisson sys- 
tem if it is a minimum of the energy functional VV at fixed 
mass. If a minimum exists, it is necessary that u = 0. 
As a result, a steady state of the Euler-Poisson system 
is nonlinearly dynamically stable if, and only if, it is a 
minimum of VV[p] = W[p, u] — 9[u] at fixed mass M . In 
conclusion, the condition of nonlinear dynamical stability 
can be written 



Min {W[p\ | M\p] = M}. 



(C4) 



For a white dwarf star at zero temperature (T = 0), ac- 
cording to Eqs. (|C2[) and (|C3[) . the functional VV reduces 
to the energy E = U + W, where U is the kinetic energy. 



2. Free energy of self-gravitating fermions 

We can also view a white dwarf star as a gas of self- 
gravitating relativistic fermions at statistical equilibrium 
in the canonical ensemble. It is thermodynamically sta- 
ble if, and only if, it is a minimum of free energy F at 
fixed mass M. The free energy is given by 





F = E - 


m J 


f {L ln L 



TS = / fe(p)dpdr 



1- 
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In 1- 



Vo 



)} 



p$dr 
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where tjq = 2/h is the maximum value of the distri- 
bution function fixed by the Pauli exclusion principle. 
For a white dwarf star at zero temperature (T = 0) , the 
free energy F reduces to the energy E = U + W . The 
condition of thermodynamical stability in the canonical 
ensemble can be written 



Min {F[f} | M[f] = M}. 



(C6) 



The critical points of free energy at fixed mass, deter- 
mined by the variational principle SF — aSAI = 0, cor- 
respond to the relativistic mean-field Fermi-Dirac distri- 
bution 



/ = 



Vo 



o /3[e(p)+nH$(r) + \ ) 



(C7) 



where A is a Lagrange multiplier (related to a) deter- 
mined by the mass M and e(p) is given by Eq. ((4|). Using 
Eqs. @, © and (|C7[) . the density and the pressure are of 
the form p = p(pH$(r) + A ) and P = P(pH$(r) + A ). 
Eliminating p_ff $(r) + Ao between these two expressions, 
we find that P = P(p) so that the gas is barotropic. The 
equation of state is parameterized by T and is fully de- 
termined by the entropic functional in Eq. (|C5[) . which 
is here the Fermi-Dirac entropy. At T = we obtain the 
explicit relations of Sec. HU Furthermore, the condition 
that /(r,p) is a function of the energy e = e(p)+pH<f>(r) 



implies that the corresponding barotropic gas is at hydro- 
static equilibrium. Indeed, using Eq. ([B]), 

VP = ~ J ^pe'(p)dp = M #V$i J f'(e)pe'(p)dp 

= M if V $i I (p ■ |^ J dp = -pHX7<S> J fdp, 



so that 



VP = -pV$. 



(C8) 



(C9) 



The relativistic mean-field Fermi-Dirac distribution 
(|C7|) is just a critical point of free energy at fixed mass 
M. To determine whether it corresponds to a true mini- 
mum of free energy, we can proceed in two steps. We first 
minimize F[f] for a fixed density profile p(r). Since the 
potential energy W[p] and the mass M[p] are entirely 
determined by the density profile, this is equivalent to 
minimizing U[f] — TS[f] at fixed p(r), where U is the 
kinetic energy. This gives a distribution 



/ 



Va 



,/3[e(p)+A(r)] ' 



(CIO) 



where A(r) is a local Lagrange multiplier determined by 
the density p(r), using p = pH J fdp. Since 5 2 (U — 
TS) > 0, the distribution (|(J10|) is a true minimum of 
F[f] at fixed p(r). Substituting the optimal distribution 
function (JCIOP in Eq. (IC5[) . we obtain a functional F[p] = 
F[f] of the density p(r). Using Eqs. ©, © and (|CT0| . 
we note that the density and the pressure are of the form 
p(r) = p(A(r)) and P(r) = P(A(r)). Eliminating A(r) 
between these expressions, we find that P = P{p) where 
the equation of state is the same as the one determined 
from the Fermi-Dirac distribution (|C7|) at equilibrium. 
Now, we can show that 



F[p] = W[p], 



(Cll) 



where W[p] is the functional defined in Sec. IC 11 The 
relation (|C11|) can be obtained by an explicit calculation 
which extends the proof given in Appendix B of [25( for 
classical particles (this will be shown in a future work; see 
also particular cases in Sec. IC 3|) but it also results from a 
straightforward argument. We note that the distribution 
function (|C10p corresponds to a condition of local ther- 
modynamical equilibrium with uniform temperature and 
zero average velocity. Thus, it locally satisfies the first 
law of thermodynamic d(u/p) = —Pd(l/p) + Td(s/p) 
with dT = 0. Integrating this relation like in Sec. IC 1[ 
we find that U[p] - TS[p] = Wi[p] where U[p] = U[f] 
and S[p] =S[f]. This directly yields the identity |(CTT|) . 
We can now conclude that F[f] has a minimum /* (r, p) 
at fixed mass M if, and only if, F[p] — W[p] has a min- 
imum p*(r) at fixed mass M. In that case, /*(r,p) is 
given by Eq. (IC10|) where A*(r) is determined by p*(r), 
writing p* = pH J f*dp. Therefore, a system is ther- 
modynamically stable in the canonical ensemble if, and 



19 



only if, the corresponding barotropic gas with the same 
equilibrium distribution is nonlinearly dynamically sta- 
ble with respect to the barotropic Euler-Poisson system. 
Said differently, a system that minimizes the functional 
(jClip is (i) thcrmodynamically stable in the canonical 
ensemble and (ii) nonlinearly dynamically stable with re- 
spect to the barotropic Euler-Poisson system. This result 
applies to white dwarf stars at arbitrary temperature. 



3. Application to white dwarf stars at T = 

Although the above results are general, it may be use- 
ful to explicitly compute the functionals (|Cip and (|C5|) 
for white dwarf stars at T = and check the relation 
(|C11[) . If we view a white dwarf star as a barotropic gas 
described by an equation of state P(p), its energy W can 
be written 



with 



1 V - / pT{p)dr + i / fxhlr. 



with 



t(p) 



p P(p') r] , 



(C12) 



(C13) 



A white dwarf star is nonlinearly dynamically stable with 
respect to the Euler-Poisson system if it is a minimum 
of W at fixed mass M. At T = 0, the equation of state 
is given by Eq. ([5]). In the classical limit (polytrope 
Tl 3/2 = D/2) we have 



W = Ml f p(D+2)/D dj . + 1 f p ^ 



(C14) 



and in the ultra-relativisitic limit (polytrope n' 3 = D) we 
have 

W = DK 2 I P (D+1)/D dr + i J p<$>dv. (C15) 

In the general case, using Eq. ((SJ), the function (|C13|) can 
be written 



rw ^r^. 



B / x 



r'D+1 



(C16) 



Integrating by parts and using Eqs. (9| and (fT0| . we find 
after straightforward calculations that 



pH 



mc 



np) = VTT^-±J o j^rj.dt. 



(C17) 



Alternatively, we can view a white dwarf star as a gas 
of self-gravitating relativistic fermions at statistical equi- 
librium in the canonical ensemble. At T = 0, its free 
energy F = E — TS coincides with its energy E. Using 
Eq. (|C5j) . it is given by 

F = E = J JL- K { p )dv+ l -J p$dr, (C18) 



fiH 



«(p) = / fe(p)dp. 



(C19) 



A white dwarf star at T = is thermodynamically stable 
if it is a minimum of the energy E at fixed mass M. In 
the classical limit, e = p 2 /2m and P = -^ J(f/m)p 2 dp 
so pn/fiH = (D/2)P. Therefore, the free energy can be 
written 



F = — I Pdv 



1 



p$dr. 



(C20) 



Using Eq. ((25|), this is equivalent to Eq. (|C14|) . In the 
ultra-relativistic limit, e = pc and P — t> J fpedp so 
pn/ pH = DP. Therefore, the free energy can be written 



F = D I Pdv + - J p$dr. 



(C21) 



Using Eq. (|3"9"|) . this is equivalent to Eq. (|C15[) . In the 
general case, we have 



n(p) 



2S D 
nh D 



fPO 

/ e{p)p D - 1 dp. 
Jo 



(C22) 



Using Eq. (fj|, it can be put in the form 

^4 = 4 f Vl + *t D -Ht. (C23) 

Now, it is straightforward to check that the two expres- 
sions in the r.h.s. of Eqs. (|C17j) and (|C23j) are equal so 
that r(p) = k(p)/ p,H implying the relation (jClljl . 

Finally, if we consider a classical isothermal self- 
gravitating gas at temperature T with an equation of 
state P — phsT/m, its free energy (1C5J) can be written 



F = —Nk B T + k B T [ — In — dr + - I p$dr. (C24) 
2 mm 2 ' 



Comparing with Eq. (|B7[) , we find that relation (jClip is 
indeed satisfied. Therefore, a classical isothermal self- 
gravitating gas that minimizes the functional (|B7[) or 
(JC24I) at fixed mass is (i) thermodynamically stable in 
the canonical ensemble and (ii) nonlinearly dynamically 
stable with respect to the barotropic Euler-Poisson sys- 
tem. We had already made this observation in 34j. In 
[2a ] , this result was extended to an arbitrary form of en- 
tropic functional, including the Fermi-Dirac entropy for 
non-relativistic systems. The present paper shows that, 
due to relation (|C11[) . the equivalence between nonlin- 
ear dynamical stability with respect to the barotropic 
Euler-Poisson system and thermodynamical stability in 
the canonical ensemble is general. 



APPENDIX D: NON VIABILITY OF A D > 4 
UNIVERSE 

If we consider, in a Z?-dimensional universe, a Hamil- 
tonian system of self-gravitating classical point masses 
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whose dynamics is described by the Newton equations 

Gm (rp - r a ) 



r» = 



E 



|r/3 



\D 



fee, 

the scalar Virial theorem reads [51] : 



(Dl) 



(D2) 



where K is the kinetic energy and Wu the trace of the 
potential energy tensor for the TV-body system [51|. For 
D 7^ 2, using Eq. (|B32[) and introducing the total energy 
E = K + W, which is a constant of the motion for an 



isolated system, we have [511 ]: 



We note that the dimension D = 4 is critical. In that 
case, I = 4E which yields after integration / = 2Et 2 + 
Cit + C-2- For E > 0, I — > +oo indicating that the system 
evaporates. For E < 0, I goes to zero in a finite time, 
indicating that the system forms a Dirac peak ("black 
hole") in a finite time. More generally, for D > 4, since 
(D - 4)1T < 0, we have I < 2Et 2 + C x t + C 2 so that 
the system forms a Dirac peak in a finite time if E < 0. 
Therefore, self-gravitating systems with E < are not 
stable in a space of dimension D > 4. The study of 
the present paper indicates that this observation remains 
true if quantum (Pauli exclusion principle for fermions) 
and relativistic effects are taken into account. In this 
sense, a universe with D > 4 is not viable. 



-I = 2K +{D-2)W = 2E+{D-4)W. (D3) 
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